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Abstract: We propose a new technique for measuring the polarization of hadronically 
decaying boosted top quarks. In particular, we apply a subjet-based technique to events 
where the decay products of the top are clustered within a single jet. The technique requires 
neither 6-tagging nor VF-reconstruction, and does not rely on assumptions about either the 
top production mechanism or the sources of missing energy in the event. We include 
results for various new physics scenarios made with different Monte Carlo generators to 
demonstrate the robustness of the technique. 
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1. Introduction 

The top quark, with its large coupling to the Higgs sector, plays an important role in models 
of physics beyond the Standard Model. Indeed, many such models posit the existence of 
top partner states (e.g. the stop squark of SUSY [1] and the T' of little Higgs models [2]) or 
otherwise couple the top to new physics in a special way (as with KK-gluons [3]). Measuring 
the couplings of the top to new states is therefore essential in distinguishing the correct 
model of physics beyond the Standard Model. 

One especially interesting aspect of these couplings is their chirality: whether or not 
they distinguish left- from right-handed tops. Fortunately, the large mass of the top, 
which makes its study so interesting for electroweak physics, makes it possible to imagine 
measuring the chiral couplings of the top directly. Unlike the other quarks, the top decays 
before hadronization, so information about its spin is transferred to the distributions of its 
decay products [4]. On the other hand, the large mass of the top also means that in order 
for the chiral couplings of tops to new physics to translate into observable top polarization 
signals, the tops must be significantly boosted, as chirality only becomes equivalent to 
helicity in the massless limit. Boosted tops are therefore a natural and interesting place to 
look for polarization signals. 

Conventional methods for measuring the polarization of non-boosted tops begin by 
reconstructing the top rest frame and considering the angular distributions of its decay 
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products in that frame, and often focus on the semi-leptonic decay mode, which can be 
fuUy reconstructed if the only missing energy in the event comes from the neutrino. Such 
techniques have been extended to events where the hadronic top is boosted, but the lepton 
from the leptonic top decay is still isolated. This isolated lepton can then be used to 
measure the polarization of its parent top, either by reconstructing the tt system [5] or 
through the shape of the lepton px spectrum [6] . 

When the top quark is highly boosted, however, requiring an isolated lepton begins to 
require a significant acceptance price. Moreover, while the large spin analyzing power of 
the lepton in standard model top decay makes it particularly useful for top polarization 
studies, it is also desirable to develop techniques which can measure polarization in boosted 
tops without the need for an isolated lepton. Being able to study polarization in boosted 
hadronic tops increases acceptance, and has the additional feature of flexibility: unlike 
leptonic tops, hadronic tops are fully reconstructable in events with multiple sources of 
missing energy. For highly boosted tops, the finite angular resolution of the detector 
makes complete reconstruction of the system difficult, and angular distributions in the top 
rest frame are no longer optimal observables. 

Here we present a technique to measure the polarization of a boosted top in its hadronic 
decay mode using the energy fraction distribution of a particular subjet. This new method 
does not require high-p^ ^-tagging, which is known to be challenging. We also do not 
require W reconstruction inside the top jet. Again, as we are considering hadronic tops, this 
technique measures top polarization using information from the top jet alone, independent 
of other objects in the event. It does not involve the reconstruction of top rest frame, or 
rely upon the measurement of missing momentum. 

While identification of boosted hadronic tops above the QCD background is challeng- 
ing, many promising approaches have been developed [7, 5, 8, 6, 9, 10, 11, 12, 13, 14, 15, 16, 
17, 18, 19]. In this article we will assume that the boosted top candidates can be identified 
through one of these methods. 

We will begin by motivating our choice of a subjet-based technique for studying the 
substructure of a top-jet. Then we will propose an algorithm useful for measuring the top 
polarization and discuss its interpretation. Finally, we will demonstrate the robustness of 
the algorithm by testing it in different physics scenarios with data from different parton 
shower models. 

2. Looking Inside a Top Jet 

Here we will discuss the different techniques used to study boosted hadronic tops. This 
will give us a chance to motivate our use of subjets while outlining other possibilities. 

In the past, two distinct approaches have been taken to analyze top jets. One approach 
uses jet shape variables [15, 16, 13] to define a function on the constituents of a top jet (in 
practice, the constituents will be calorimeter cells), treating each constituent independently. 
The other approach [13, 14, 18, 19] defines a function on the subjets formed by reclustering 
the constituents of a larger jet. Functions then depend upon the constituent four-momenta 
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only through the total four- momentum of the sub jet they are clustered into, rather than 
upon each constituent four-momentum independently. 

Each approach has both advantages and disadvantages. Subjets can reduce the effects 
of soft contamination^ by summing together constituents so that softer particles have a 
proportionally smaller influence. However, care must be taken because some quantities one 
can form from subjets, such as invariant mass, can be extremely sensitive to calorimeter 
spacing and out-of-cone radiation. Fortunately, as long as one avoids these troublesome 
quantities a subjet-based analysis can be made fairly robust. For our algorithm below, 
we will only rely upon the relative hardness and separation of the subjets, both quantities 
which are fairly insensitive to additional soft radiation and detector effects. 

Jet shape variables, because they treat each jet constituent independently, are more 
amenable to higher order calculations than variables defined with subjets. However, these 
variables can become very sensitive to the effects of contamination. As an example, consider 
the planar flow iet shape of [15], which is equivalent (up to an overall constant) to det S"*- 
defined in [13]. The planar flow of a jet is defined as 

(Ai + 

where Ai^2 are the two eigenvalues of the matrix 

/^' = E-^ — - (2-2) 

where Wi is the energy and Pi^k the fcth transverse momentum component of the ith jet con- 
stituent. This quantity essentially decomposes the jet's radiation into two moments \i^2^ 
similar to moments of inertia, so that if the jet is symmetric about its center then Pf ~ 1. 
Planar flow is useful in top-quark studies because top jets are relatively symmetric about 
their center (corresponding to higher values of planar flow), while QCD events are domi- 
nated by a single emission (corresponding to a lower value of planar flow). Unfortunately, 
planar flow weights each constituent according to its transverse momenta relative to the jet 
axis, so that as the radius of a jet is varied soft radiation towards the edge of a jet begins 
to dominate and all events are skewed toward higher Pf^. To demonstrate this sensitivity 
and how it can be reduced through the use of subjets, we have included Fig. 1, showing the 
calculation of planar flow at matrix element level, after showering, and after reclustering 
using subjets. The subjets are formed using the procedure described in Section 4 using 
R = 0.2 cones. Here one can see the large corrections to the matrix element results that 
are attributable to soft radiation. Of course, one can mitigate this effect by using smaller 
cones (the authors of [15] used R = 0.4), as the amount of diffuse soft radiation clustered 
into the top jet goes roughly as R? , yet even in this regime the effect of soft contamination 
can still be significant, especially near Pf ~ 0. 



^Contamination, radiation clustered within the top jet that did not arise from the top decay, can be the 
result of initial state radiation, multiple interactions, or wide angle emissions from other parts of the event. 

^While soft effects can skew Pf toward higher values, it remains difficult for QCD jets to be skewed all 
the way toward Pf ~ 1. Thus, while the distributions of Pf are sensitive to soft effects, they can still be 
used as effective top/QCD discriminants. 
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Figure 1: Comparison of planar flow for left-handed tops (left figure) and right-handed tops (right 
figure). For each set of samples we compute the planar flow using the three partonic decay products 
of the top (black, solid), the constituents of the showered jet using a small R — 0.4 cone (blue, 
dotted), the constituents of the showered jet using a larger R = 1.0 cone (magenta, dot-dashed), 
and the three i?sub = 0.2 subjets formed from the R = 0.4 top-jet constituents (red, dashed). These 
events are taken from the decay of a 3 TeV Z' into two tops, clustered using the anti-fcj^ algorithm, 
where we have required the top jet's mass satisfy 140 GeV < mj < 210 GeV. Note that the peak of 
the red/black distributions near Pf ^ can be ascribed to configurations where two of the partons 
become coUinear, or where one parton becomes particularly soft. 



To be sure, jet shape variables (including planar flow) are very useful and will likely 
play a role in boosted top chirality measurements. A simple counting exercise shows that 
after requiring the reconstruction of the W mass, top four-vector, and allowing for axial 
symmetry, there are still three remaining degrees of freedom encoded in the matrix element 
that can be mapped out by jet shape variables. However, to simplify the discussion and 
avoid complicated issues of contamination and higher-order corrections we will use the rest 
of the paper to focus on quantities calculated using subjets. 



3. Top Polarimetry With Subjets 

We will now explore methods for using subjets to measure the polarization of a coUimated 
hadronic top. In what follows, we will assume we are working with jets tagged as tops, 
as discussed in [13, 15, 16, 14, 20], and subsequently decomposed into three subjets (a 
prescription for such a decomposition is given later). 

3.1 Choosing a Polarimeter 

One observable sensitive to the polarization of the top is the distribution of energy among 
the its three decay products in the lab frame. In the collinear limit, the lab-frame energy 
fraction of the ith subjet, Zi = Ei/Etop, depends only on the energy and angular distri- 
butions in the top rest frame, and can serve as a robust variable to measure polarization. 
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While energy fractions are not Lorentz invariant for finite top mass, and in particular are 
not invariant under longitudinal boosts, frame dependence enters only at order mt/Et, 
and therefore, for highly boosted top quarks, energy fraction variables become fixed, sta- 
ble quantities ^. The question then becomes how to select the subjet to be used as a 
polarimeter. 

The most obvious candidate for the job is the 6-jet [21, 22], identified either directly 
through 6-tagging or indirectly by first finding the light quarks from the W. However, the 
identification of the b and W poses some experimental difficulties. Even in isolation, the 
efficiency of 6-tagging drops by a factor of 2-3 at high px while light quark rejection is 
degraded by roughly a factor of 3 [6, 10, 23, 24, 25]. When the 6-jet is situated within a 
collimated top jet, the additional tracks from the neighboring light quark subjets present 
added complications for 6-tagging algorithms. 

Another possible method of identifying the 6-jet is to do so indirectly, by finding the 
W. One approach to identifying the W is to look for two jets with an invariant mass 
within the W mass window. However, the subjet invariant mass distributions are distorted 
both by contamination from soft radiation and by imperfect subjet reconstruction, as well 
as by the finite size of the calorimeter. The invariant mass mij of two nearby subjets is 
approximately proportional to their separation in R, and for subjets whose centers are 
separated by ARij < 0.5, the uncertainty associated with the calorimeter granularity 
6R ~ 0.1 can be significant. Distinguishing the correct W subjet pair from amongst the 
three choices, all of which are typically within a factor of two of each other, then becomes 
difficult. 

Another possible strategy to identify the 6-jet is to look for hard splittings within the 
top jet. As discussed in [13], the energy sharing of a parton branching A — t- BC, 

z{A ^ BC) = mm{EB, Ec)/Ea, (3.1) 

discriminates between hard splittings from decays, z ~ 0.5, and soft splittings more char- 
acteristic of QCD, 2 ~ 0. If the W decay products were well-separated from the 6-jet, one 
could identify the b by unwinding the clustering of the top jet until there were two subjets 
and tagging the b as the one with smaller z (so the W subjet would be the one with a 
harder splitting). Unfortunately, because the W has a mass on the same order as that of 
the top, the W decay products are not well-separated from the 6, so upon unwinding the 
top jet by one step one often finds that the &-jet has been clustered with a lighter jet from 
the W decay. 

We propose here an alternate subjet selection algorithm, based on kx distances be- 
tween subjets, which does not require either b or W identification. While the algorithm 
is conceptually less straightforward than those based on attempting to identify specific 
partons within the top jet, it yields a distinct separation between chiralities and is ro- 
bust under showering and detector effects. Consider the distance measure between two 
four- momenta i and j, 

dij = mm{pli,p^j)Rij, (3.2) 

^Depending on the boost of the top quark, it might also be desirable to consider subjet pr fractions, as 
the mt/Et corrections to the collinear limit differ for energy and pt fractions. 
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Figure 3: Energy fraction of the parton selected by our algorithm, broken down by parton identity. 
The events shown here correspond to tops produced from a 3 TeV resonance. 



;)^. Of the three dij one can form from the three top 



where Rfj = {rji - Vj)'^ + i4>i ■ 
subjets, consider the smallest. Our top polarimeter is the energy fraction zk of the harder 
jet jx in the minimum distance pair. We plot the distribution of this variable at parton 
level for different chiralities in Fig. 2. The variable shows a clear distinction between right- 
and left-handed top quarks, with right-handed tops peaked at smaller values of zk, and 
left-handed tops preferring larger values of zk ■ 
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Figure 2: Energy fraction zk of the parton se- 
lected by the fcy-based algorithm for different top 
polarizations. The events shown here correspond to 
tops produced from a 3 TeV resonance. 



3.2 Operation of the algorithm 

The success of the jet Jk selected by this 
algorithm polarimeter depends on 
multiple aspects of the angular and en- 
ergy distributions of daughter partons in 
polarized top decay, which for reference 
are reviewed in the Appendix. In order 
to explain the success of our polarime- 
ter, we first consider how the algorithm 
functions at parton level. The identi- 
ties of the partons picked out by the al- 
gorithm differ between right- and left- 
handed tops. In Fig. 3 we break down 
the contributions to the variable zk by 
parton identity. 

The anti-down quark is maximally 
correlated with the top spin, and thus 
for left-handed tops the d tends to be 
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soft. For left-handed tops the minimum- 
pair therefore tends to involve the d, and in such pairs the other parton (6 or u) is the 
harder of the two. The algorithm therefore picks out first b quarks, which take a larger 
fraction of the top energy, and secondarily u's, with d quarks a distant third. 

For right-handed tops, where the top energy is shared more equitably among the 
daughter partons, angular correlations play a more central role. The d-quark is now both 
more central and harder than predicted by pure phase space (due, again, to its maximal 
correlation with the top spin). Therefore in order to reconstruct the necessary invariant 
masses, the AR separation between the d and the u and b quarks must be smaller than 
for pure phase space, and the minimum /cj- pair then tends to involve the d. For right- 
handed tops, the algorithm thus dominantly selects the d-quark, as can be seen in Fig. 4. 
While the d is preferentially emitted along the top direction of motion, its energy fraction 
distribution nonetheless falls off at high energies, as the lab-frame d-quark energy fraction 
depends on the energy of the d-quark in the top rest frame as well as the angle of emission. 
The contribution of the 6-quark to the variable zk comes mostly from hard 5's recoiling 
against soft transverse VF's. 

At high parton energy fraction zk, the algorithm dominantly selects the hardest par- 
ton: b and u for left-handed tops, b and to a lesser extent u, d for right-handed tops. At 
intermediate energy fractions, the origin of the u and d partons from a common W comes 
to dominate. The distance between the u and the d is bounded from above, as the u 
and the d must reconstruct the W. In events without hierarchical energy distributions, the 
minimum k^ distance thus tends to be between the decay products of the W. Therefore, 
at intermediate energy fractions, the parton selected by the algorithm is predominantly 
the d (for right-handed tops) or the u (for left-handed tops). This can be seen already in 
Fig. 3, and is further demonstrated in Fig. 4. 

Finally, we note that all of these arguments are based upon the assumption that one can 
go from the collider coordinate system to one oriented around the top direction of motion 
without significant effects. This assumption does not hold exactly, because the detector 
geometry is not invariant under rotations around the axis defined by the top direction of 
motion, and because the kj- algorithm used to select the subjet jx makes reference to 
the collider coordinate system through the definition of transverse momentum. Therefore, 
events which differ from each other only by a rotation around the top axis of motion appear 
differently both in the detector and in the subjet selection algorithm. Interference terms 
between right- and left-handed tops generically then do not completely cancel. However, 
as the magnitude of the interference contribution is determined by the components of the 
parton momenta transverse to the top momentum, these effects are of order mt/Et, a 
subleading effect for large boosts. 

3.3 Implementation 

To implement this algorithm in practice one must have a technique for finding three subjets 
within the top jet. The exact procedure one uses to identify the subjets is not important, 
but for concreteness we detail the method used in our study. One advantage of our subjet- 
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Figure 4: Energy fraction of the b quark, broken down according to which two partons in the 
event belong to the pair with minimum kx for left-handed tops (left side) and right-handed tops 
(right side). Note that the contents of the plots are stacked. At high energy fractions, the W and 
its decay products are soft, and the minimum fc-r pair tends to involve one of the W decay products 
together with the b. At intermediate energy fractions, the minimum kx pair tends to be the W 
decay products. At small energy fractions, the b begins to appear as the softer of the two partons 
in the minimum kx pair. The effect is more pronounced for right-handed tops, which have a less 
hierarchical distribution of energy among the three daughter partons. The events shown correspond 
to tops produced from a 3 TeV resonance. 

finding technique is that it is easy to implement within the FastJet [26] framework, already 
used by many studies for top tagging. 
The procedure is as follows: 

• Cluster the event with a reasonably sized cone {AR > 0.7) and select a top candidate. 

• Take all the cells clustered into the top candidate and recluster them using a smaller 
cone {AR ^ 0.2). 

• Demand that there are at least three subjets, each with a substantial amount of the 
jet's energy ^ 1 — 2%. If there are not, split the harder subjet by unwinding [13] it 
one step using the kx algorithm, and use the two resulting daughters along with the 
second hardest subjet from before. 

• Now use the four-momenta of the three subjets to find the pair with the smallest kx 
distance measure, and compute the ratio of the energy of the more energetic jet in 
this minimum- A;7^ pair to the energy of the entire top jet. 

Results using this procedure are shown in the next section. 
4. Examples 

We will now apply the subjet-based technique developed in the previous section to some 
realistic examples. Our goal is to show that the technique works for fully showered events 
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Figure 5: Energy fraction of the selected jet for results from (left to right) Herwig++, Pythia-6 
(Q^) and Pythia-6 (pt)- For each plot the solid red and dotted black lines come from right- and 
left-handed tops, respectively. 

clustered with finite calorimeter cells using a variety of parton shower and hadronization 
algorithms. It is important to note that we do not consider the shape of background QCD 
distributions, nor do we consider any shaping effects that might arise from the effects of 
top tagging. A more complete experimental study would include these effects, but due to 
the high discriminating power of top tagging algorithms (not to mention other aspects of 
the event that could be used to remove background) and their relatively high efficiency, we 
do not expect these effects to be significant. 

In what follows, our analysis is performed on events generated at matrix element level 
using MadGraph 4.4. 17 [27] for physics at the LHC scale (14 TeV). Subsequent showering 
and hadronization is performed using Pythia 6.4.21 [28] and Herwig++ 2.3.2 [29] . When 
using Pythia, we consider parton showers generated using both virtuality (labeled as Q^) 
and Pt ordered showers. Visible final state particles are grouped into 0.1 x 0.1 calorimeter 
cells before being clustered into R = 0.7 jets using the smti-kT [30] algorithm. To ensure 
that the top decayed into visible products (and that no significant radiation was lost outside 
the cone) we demand that the jet mass exceed 170 GeV. We form subjets by running the 
anti-Zsr algorithm with R = 0.2 on the constituents of the top jet, requiring that the third 
most energetic subjet carry at least 1% of the top jet energy, and splitting the hardest 
subjet if this condition is not satisfied. 

4.1 Tops from a resonance 

We begin by studying a colored octet vector G' with a chiral coupling to the top quark. 
This model was chosen for simplicity, but it captures the main features of well-motivated 
scenarios like KK-gluon production. The process under consideration is 

gg^G'^ tt (4.1) 

where mc = 3 TeV. The results are shown in Fig. 5. One can see from these distributions 
that the characteristic shapes from matrix element level are unchanged after parton show- 
ering and hadronization, demonstrating the robustness of our subjet selection technique. 
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Figure 6: Energy fraction of the selected jet for results from (left to right) Herwig++, Pythia-6 
(Q^) and Pythia-6 (pt)- For each plot the solid red and dotted black lines come from right- and 
left-handed tops, respectively. 

4.2 Tops from cascade decays 

Cascade decays of an on-shell top partner (such as a stop squark or a T') to a top plus 
missing energy are a standard signal of a broad class of well- motivated models. In the 
presence of multiple sources of missing energy, there is no longer enough information to 
solve for the rest frame of a leptonically-decaying top quark. Hadronic top quarks, which 
can be reconstructed independently of the other particles in the event, become a more 
useful source of information. 

For tops produced from an un-reconstructed parent, the observable polarization signal 
depends on the masses of the new physics particles through the relation of the unknown 
parent rest frame to the lab frame, as well as through the vertex kinematics [21]. The lack 
of information about the parent rest frame reduces the observable polarization signal for 
tops coming from a cascade decay compared to the signal from a resonance. Nonetheless, 
observable signals are still possible as long as the boost of the top from its parent is 
dominant, allowing the chiral structure of the top production vertices to be probed. 

We consider a model for production of two top partner T' particles decaying into tops 
and sources of missing energy (labeled A^). Our choice of spectrum has mx' = 2 TeV and 
m^o = 100 GeV. The results of the analysis performed on the model are seen in Fig. 6. 
As for the G' , the distributions agree well with the parton level results and have the same 
qualitative shapes regardless of the generator used. 

5. Conclusions 

We have proposed an analysis tool useful in determining the chiral structure of the top 
quark's coupling to new physics. Our method uses subjet-based techniques to probe sce- 
narios where a highly boosted top decays hadronically. This tool requires no assumptions 
to be made about the production mechanism of the top or about the origin of missing 
energy in the event, and does not rely upon 6-tagging or W reconstruction. 

By testing our method on Monte Carlo data from multiple generators using different 
new physics scenarios we have indicated its robustness against the effects of parton show- 
ering and calorimeter segmentation. A more complete analysis would study the shaping of 
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the distributions from the top tagging method using in selecting a sample, but we expect 
these effects to be small. 

Boosted hadronic tops may provide a new window to shed light on otherwise difficult 
aspects of new physics at the LHC, and will certainly provide a complementary probe 
of physics beyond the Standard Model. Variables which can analyze the polarization of 
boosted hadronic tops, such as those introduced here, will fill an important slot in the 
analysis toolkit as we try to unravel the physics behind LHC data. 
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A. Angular Distributions in Decays of Polarized Tops 

Here we collect some results on the energy and angular distributions of daughter partons 
in polarized top decay. To arrive at these results we assume the standard model top decay 
t —7- Wb with subsequent W decay W — )■ du, and work at leading order in the narrow width 
approximation for both t and W. Then, the squared matrix element for a polarized top 
can be written 

\M\^ oc {t± ■ d){b ■ u) (A.l) 

where d, b, and u denote the momentum four-vectors of the respective partons, and t± = 
Pt ± mtS contains information about the top polarization through the spin four-vector S. 
In the top rest frame, the spin four-vector takes the form = (0,s), where s is a unit 
vector defining the axis of polarization. In the narrow width approximation, and further 
taking = 0, the energy of the b quark is fixed at Ej, = {ml — m^^)/2mt in the top rest 
frame. The full differential decay rate then depends nontrivially on two quantities, which 
we can take to be the down-quark energy in the top rest frame, and the angle of the 
down quark with respect to the top spin axis, cos^^. The differential decay rate can be 
written 

1 dV 12m^ 

r^Tp^ fr= 6 % 4 2 ^d{mt-2Ed){l + Vtcos9d) (A.2) 

i dEdd cos dd rrq + 2m^^ — irriy^mf 

where Vt is the top polarization, — 1 < "P* < 1 with Vt = 1,-1 corresponding to right and 
left handed tops, respectively. 

These results can be generalized to study the angular distribution of any given daughter 
parton i in the top rest frame by writing [31] 

'l + Vt^iCosOi), (A.3) 



r d cos 9i 2 

where Kj is the spin analyzing power of the parton z, and cos^j is the angle that parton 
makes with respect to the top spin axis in the top rest frame. The values of for various 
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Parton 


h W d u j 


K 


-0.4 0.4 1.0 -0.3 0.5 



Table 1: Tree level values for the spin analyzing power k of various top daughters in top decay. 
The object j is defined in the text. 

choices of i are listed in Tab. 1. The d, which corresponds to the lepton in leptonic top 
decays, is maximally correlated with the top spin, with = 1. In addition to the partonic 
b, W, u, and d, we also consider an object j, defined to be the softer of the two light quark 
jets in the top rest frame. As the d tends to be softer than the u in the top rest frame, 
this jet is the d-jet approximately 60% of the time [32, 33]. 
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